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Abstract 

In this talk I report on the results of a complete 0{a) calculation of leptonic 
//-decay. The calculation is complete in the sense that all polarization and mass 
, effects have been included in the radiative corrections. I mostly concentrate on the 

I longitudinal polarization of the electron which considerably differs from the naive 

' value Pj = — 1 in the threshold region, both for the Born term and more so for the 

■ radiatively corrected case. I also discuss the role of the 0(a) anomalous spin-flip 
D I contribution and its description in terms of the equivalent particle approach which 

" ■ survives the rUe limit. Finally, I provide a brief account of the history of the 

. 0(q) radiative corrections to leptonic /i-decay. I trace the error done in the first 

I (wrong) 1956 calculation of Behrends, Finkelstein and Sirlin. My account of this 

T-H ' historical error differs from that recently given in a talk by Kinoshita on the occasion 

■ of the 70th birthday of Sirlin. 

m ■ 
(N ■ 

^ ' 1 Introduction 

O ' Let me first clarify what I mean by "complete" 0{a) solution of the /i-decay problem. 

r-| ■ With polarization and mass effects included, the decay /i~(T) + + Ue is de- 

0-{. scribed by five structure functions Gi,..,G5. Here I have disregarded a sixth T-odd 

OhI structure function Gq which is zero in the Standard Model. The five structure functions 

are associated with the rate (Gi), single spin effects (G2 for the muon, G3 for the electron) 
and spin-spin correlation effects (G^jG^). The complete 0(a) solution of the //"-decay 
problem consists in calculating all five structure functions at next-to-leading order (NLO) 
in differential and integrated form keeping rUe 7^ 0. 
5^ I In most previous calculations the strategy of calculating the radiative corrections to the 

five structure functions was characterized by the statement "put the electron mass to zero 
whenever possible", keeping the mass dependence only in logarithmic factors involving 
ln(me/m^). This holds true for the first correct calculation of the unpolarized spectrum 
function Gi given in y]. The corresponding mg 7^ result can be extracted with a little 
bit of labour from the papers of [21 . 

Putting the electron mass to zero (whenever possible) is really quite adequate for 
most of the electron spectrum since the ratio (me/m^) is quite small. This is no longer 
true for the threshold region where electron mass effects become important for /5-values 
smaller than f3 ^ 0.995. This will be discussed later on when we calculate the longitudinal 
polarization of the electron. Moreover, after the discovery of the r-lepton it was not so 
clear to what extent the mass of the final state lepton could be neglected since in the case 
r~ ^ H~ + Ut- + the mass ratio {171^/111^-) is not really very small. 

Historically, the program of calculating the 0{a) corrections to the five structure 
functions proceeded in several steps. Here one has to distinguish between results for the 
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differential energy distributions which are characterized by the spectrum functions Gj, and 
the fully integrated results which are characterized by the rate functions Gi = J f3xGidx. 

Concerning the spectrum functions, one can extract the rrie ^ 0{a) radiative cor- 
rections to the spectrum function Gi from |21 Ej as remarked on before. The authors of 
[H] also give results for the spectrum function G2, however, for rrie = 0. In |1] Fischer 
and Scheck calculated the rrie = radiative corrections to Gi, G2, ^3 and G4 where one 
should mention that G5 describing the transverse polarization of the electron vanishes 
in the limit rrie 0. rrie 7^ results for the polarized spectrum functions G2 and for 
G2, G3, G4, G5 can be found in and in 0, respectively. 

Results for the nif, 7^ integrated rate function Gi were first obtained by Nir in the 
context of semileptonic heavy quark decays [7]. He used a different route to obtain the 
integrated rate. He started with the differential g^-distribution calculated in P] and then 
obtained the total rate by g^-integration. Because semileptonic heavy quark decays have 
a structure similar to leptonic /^"-decays in the Standard Model his results apply also to 
leptonic /i-decays. The rrie ^ polarized rate functions G2 and 62, ^3, 64, G5 were again 
obtained by Arbuzov j3] and by us respectively. It is not quite clear what caused the 
long delay from 1958 to 2003 to complete the missing 0(a) pieces in leptonic /i-decays. 
Part of the reason is that the necessary calculations are quite involved and therefore had 
to wait for the advent of modern computers with their symbolic computation facilities. 

For us calculating the complete mg 7^ 0{a) radiative corrections to //-decay was a 
natural outgrowth of corresponding 0{as) polarization calculations that we had done for 
the semileptonic quark decays t — > b + l~^ + 1/1 and b ^ c + l~ + ui ^ which, in the Standard 
Model, have a structure similar to leptonic //-decay. These calculations had been done 
keeping the full mass dependence of the final state quark. This is important in particular 
for the case b c + l~ + ui since, in this case, the mass of the final state charm quark 
mass can certainly not be neglected. The expertise in handling the full mass dependence 
in semileptonic heavy quark decays was then exported to leptonic /i-decays. 

2 Angular decay distribution 

The angular decay distribution of the semileptonic decay of a polarized muon into a 
polarized electron is given by [H] 

-- (3x Vq{Gi + G2P cos Op + G3 cos e + G4-P cos Op cos 



dx dcosOp 

+ GsPsin^psin^cosx + Gg-Psin^psin^sinx) , (1) 

where is the polar angle between the electron and the polarization vector P of the 
/i~ in the /t~ rest system , and 9 and x are the polar and azimuthal angles describing 
the orientation of the spin quantization axis of the electron The Born term rate 
for vanishing electron mass mg = is given by Fq = Gpm^/1927r^. As usual I have 
defined a scaled electron energy x = 2Ee/m^. Further, the velocity of the electron is 



determined by /3 = \'Pf.\/Ee = v 1 ~ {'iy'^/x'^) where y = me/m^. As mentioned before. 
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Gi is the unpolarized spectrum function, G2 and G3 are single spin polarized spectrum 
functions referring to the spins of the /i~ and e~, resp., and G4, G5 and Gq describe 
spin-spin correlations between the spin vectors of the muon and electron. Gq represents 
a so-called T-odd observable which is identically zero in the Standard Model [H]. As 
mentioned before, the structure function G^ is proportional to the electron mass and 
therefore vanishes for rrie —>■ 0. 

3 Longitudinal polarization of the electron 

Here I concentrate on the longitudinal polarization of the electron which, for an unpolar- 
ized muon, is given by 

pi(x) (2) 
For the Born term contribution one obtains 

PHx) - 3 ^(3-2x + y^) 




-2 -1 
10 10 1 
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Figure 1: Longitudinal polarization of the electron in leptonic muon decays at LO (full 
line) and NLO (dashed line) as a function of the scaled energy x = 2Ee/m^ The NLO 
(dotted line) curve corresponds to keeping rrie 7^ in the Born contribution and mg = 
in the 0{a) contribution. 

The LO curve starts to deviate from the naive result = —1 at around x = 0.1 
(/3 = 0.995) where mass effects start setting in. It is curious to note that the LO curve 
is very well described by the functional behaviour Pg = —f3. The correction to the 
approximate result Pj = — /5 is of 0(1%) or less such that the correct and approximate 
LO curves are not discernible at the scale of Fig. ^ I make mention of this fact since in 
some recent text books the result Pj = —P has been claimed to be an exact result for 
left-chiral heavy fermions ^Hl- Eq-© shows that Pg = —(3 cannot be an exact result. 
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Numerically, however, Pj, = —(3 appears to be a very good approximation, at least in this 
application. 

The radiative corrections to the polarization of the electron become quite pronounced 
starting at x = 0.1 {(3 = 0.995). We also show an approximate NLO curve where rrie 7^ 
for the Born term and me = for the 0{as) contribution as suggested by 0]. The 
approximate NLO curve does not have the correct threshold behaviour ~ —jS near 
threshold at x = 2?/ ~ 0.00967 resulting from the wrong threshold behaviour of the 
radiative corrections which starts showing up at x = 0.02 (/5 = 0.875). 



4 Calculation of the tree graph contribution 

The NLO charge-side tensor describing the /x^ — > e~ + 7 tree-graph transition in charge 
retention form is given by [Hj 



7— spin 



.k-pe\ k-pe k-p^ 



I ^ k-p, + ml Pe-PA.ja0 j0a j .0. 

+ k-p,[ k.p, k-p^rp^^'^p^ ^^^^ ^ 

+ ^k.pe){k.p,) + " ""'^"'^ " (k-PeKk-p,) + *^ ~ ^"'^ 



2 / ^2 ^2 



A* I 2,pe-p^ \/-a-0 , -0-a - - a0 

+ n ^ ~ 7"; ^77^ ^PePii + PeP^ - Pe-Pt.9 ^ 



2 \{k-p^f (k-pe)'^ {k-pe){k-p^) 

The momentum of the radiated photon is denoted by A;. I have ommitted an a ^ (3 
antisymmetric piece in Eq. (@)) for the reason that it does not contribute to the differential 
electron spectrum. 

Eq.(0]) is written in a very compact way. First, polarization effects are included by 
making use of the very compact notation 

— ck cx. ct —cx. ol <y. / r \ 

P^,=P^^- , Pe = Pe - ^eSe , (5) 

where and are the polarization four-vectors of the /i~ and e~. Second, in the last 
line of I have isolated the infrared singular piece of the charge-side tensor which is 
given by the usual soft photon factor multiplying the Born term contribution. Technically 
this is done by writing 

_ f(j{a)a0 _ c(")"/5(soft photon)) + C(")"^(soft photon) . (6) 



The remaining part of the charge-side tensor in (@)) is referred to as the hard photon 
contribution. It is infrared finite and can thus be integrated without a regulator photon 
mass. The integration of the infrared singular piece with a regulator photon mass will be 
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discussed in the next section including a discussion of the errors that had been made in 
the first evaluation of the infrared contribution [2]. 



5 Historical note on soft photon regularization 

Let me now turn to the calculation of the soft photon contribution. The soft photon 
(s.ph.) transition matrix element reads 

The integration over the relevant phase space is determined by 

/ = — / dz \M,,ph.\\ (8) 

47r j-i JO fco 

where z = cos 6 is the cosine of the angle between the electron and the photon, and 
kmax{z) is the maximal photon momentum kmax{z) = m^(l + — x)/{2 — x(l — Pz)) . 
The squared soft photon matrix element in (jSI) is given by 

iM,,,,,p = y: M:^kM*/pk.<w^pw ■ (9) 

A=±l 

The spin summed product of polarization vectors in can be replaced by the metric 
tensor as follows 

^ e-(A)e^(A) = diag(0, 1, 1, 0) diag(-l, 1, 1, 1) = -g^^ (10) 

A=±l 

This corresponds to using the Feynman gauge which was also used to calculate Eq.(jH). 
The longitudinal and scalar pieces can be added to the spin sum in ()10|) since their 
contributions in Q cancel due to gauge invariance. This can be explicitly checked by 
evaluating the longitudinal and scalar pieces in the rest system of the fi~ with the z- 
axis along the photon direction, where = (fco;0,0,fc), = (m^; 0,0,0) and pe = 
{Ef,; \pe\ sin6',0, \pe\ cos 9). For the squared soft photon matrix element the replacement 
in (Uni) leads to 

It is instructive to split the squared soft photon matrix element into its transverse and 
longitudinal/scalar part. One obtains 

\M.,,f = lAwii + iM.,.ia = nj^^^ + - ■ (12) 

The longitudinal and scalar contributions in \Ms,ph.\c proportional to k^ and k"^. For 
on-shell photons with ko = k the longitudinal and scalar contributions cancel as asserted 
before. 
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After doing the azimuthal integration the soft photon factor (jSJ is given by 



I = ^ dz dk-^-^ — -— r . (13) 

Similar to Eq. ()12|l it is instructive to spht the integral into a transverse contribution It 
and a longitudinal/scalar contribution Iq. One obtains 

(14) 

At this point we introduce a (small) photon regulator mass through k\ — k?' = m^. 
The regulator mass is used to regularize the infrared singularity present in Eqs. fll3|l 
and (fT^ . One can still use the metric tensor in (fTUI) since the scalar piece /c"/c^/m^ in the 
spin 1 propagator gives zero when contracted with the square of the soft photon matrix 
element M^^y^ ^s^h ■ Technically, the integration is best done by changing to the variable 
t via k = m^{t^ — l)/2t in order to get rid of the square root factor in the photon energy 



kQ = ^A;2 + m2. MATHEMATICA will do the rest for you. 

For reasons of brevity I only list the y ^ results of the two integrations in ()14|1 . One 

has (A = m^/m^) 

A ^ ^ ' ^ ' 12 ' ^ X 



2 In i-^ ( In(-) - l) + U^{x) - + + {I ~ -) ln(l - x) - 1 , (15) 



and 



Ic = In(-) + 1 - X • (16) 

y 4 



For the sum of the two contributions one obtains 

It + Ic=( In(-) - l) (2 In - ln(-)) + U,{x) - ^ + (1 - 1) ln(l - x) . (17) 

Behrends, Finkelstein and Sirlin in their 1956 paper j,2| set ko = k in the integrand 
of the transverse contribution of ()14j) . Of course, setting ko = k in the total contribution 
fll3|l gives the same result in agreement with the arguments presented before. They then 
obtain (using their notation) 



1 Jo 'ko{l-pzy 2 7-1 il-(3zy 

= ln^i^(ln(-)-l)+L22(x) + (l--)ln(l-x)-l. (18) 
2A ^ y ^ X 

It is obvious that the result ()18p differs from the true transverse contribution ()15|). Ap- 
parently, Behrends, Finkelstein and Sirlin committed two mistakes in their 1956 paper 
First, they considered only the transverse contribution of the massive photon instead 
of including also its longitudinal part. Second, they did not correctly calculate the trans- 
verse contribution by erraneously setting k^ = k in the matrix element. Kinoshita in a 
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recent talk given on the occasion of the 70th birthday of Sirhn identifies only the first 
mistake [TT] . 

The two mistakes were corrected in the 1959 paper of Kinoshita and Sirlin by 
adding a compensation term C (their notation) to their 1956 result, where 



C = It + Ic-V=( In(-) - l) (2 In 2 - ln(-)) + 1 - ^ • (19) 

Why so much ado about a mistake done many years ago? Well, the mistake is in 
some sense historical since the 1956 result violates the so-called Kinoshita-Lee-Nauenberg 
theorem formulated a few years later on ^21 ^] • The Kinoshita-Lee-Nauenberg theorem 
states that the integrated rate should not contain any logarithmic dependence on the 
electron mass. However, when integrating the 1956 result one finds 

r(i956) = ro(^l + - vr^ + 4 [ In^y + (3 + 21n2) Iny + 2 + 41n2 + ^ ] ) ) , (20) 

which clearly violates the Kinoshita-Lee-Nauenberg theorem. The correct result is ob- 
tained by setting the square bracket in (j^n|l to zero which shows that the correct result 
does in fact satisfy the Kinoshita-Lee-Nauenberg theorem. In fact, the observation that 
there is no logarithmic mass dependence in the leptonic /i decay rate was very likely a 
progenitor of the celebrated Kinoshita-Lee-Nauenberg theorem. 



6 The anomalous spin-flip contribution 

In this section we concentrate on one interesting aspect of the /x-decay problem, namely 
on the so-called anomalous spin-flip contribution which, even in the chiral limit, flips the 
helicity of the final-state lepton at NLO. 

Collinear photon emission from a massless fermion line can flip the helicity of the 
massless fermion contrary to naive expectation. This has been discussed in a variety of 
physical contexts (see e.g. references in [6 ). This is a "me/me" effect where the me in 
the numerator is a spin flip factor and the nie in the denominator arises from the collinear 
configuration. In the limit mg the helicity flip contribution survives whereas it is not 
seen in massless QED. 

We shall discuss this phenomenon in the context of the left-chiral /x — e transition. 
At the Born term level an electron emerging from a weak {V — A) vertex is purely left- 
handed in the limit rrte = 0. Naively, one would expect this to be true also at 0{a) or 
at any order in a because in massless QED photon emission from the electron is helicity 
conserving. 

Let us take a closer look at the anomalous helicity flip contribution in leptonic — > e 
decays by considering the unnormalized density matrix element p++ of the final state 
electron which is obtained by setting cos^ = 1 in (Q) (remember that vanishes for 
me — >■ 0, and Gq = in the Standard Model). One has 

-^=(3xro{Gi + G3) . (21) 
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Contrary to naive expectations one finds non-vanishing right-handed (++) contri- 
butions which survive the mg hmit when one takes the mg hmit of the NLO 
contributions to (PT|) In fact, one finds 

5;;^.iir„((l-.f(5-2.)). (22) 

The result is rather simple. In particular, it does not contain any logarithms or 
dilogarithms. The simplicity of the right-handed contribution becomes manifest in the 
equivalent particle description of yU-decay where, in the peaking approximation, /i-decay 
is described by the two-stage process /i~ e~ followed by the branching process e" 
e~ + 7 characterized by universal splitting functions Dnf/hf{z) Pl]. The symbols nf and 
hf stand for a helicity non-flip and helicity flip of the helicity of the electron. In the 
splitting process z is the fractional energy of the emitted photon. The off-shell electron 
in the propagator is replaced by an equivalent on-shell electron in the intermediate state. 
Since the helicity flip contribution arises entirely from the collinear configuration it can 
be calculated in its entirety using the equivalent particle description. 

The helicity flip splitting function is given by Dhf{z) = az/{2Tc), where z = ko/E' = 
(E' — E)/E' = 1 — x/x', and where ko is the energy of the emitted photon. E' and E 
denote the energies of the initial and final electron in the splitting process. The helicity flip 
splitting function has to be folded with the appropriate rrie = Born term contribution. 
The lower limit of the folding integration is determined by the soft photon point where 
E' = E. The upper limit is determined by the maximal energy of the initial electron 
E' = One obtains 



dx Jo Jo dx 



dx' / dz — Dhf{z)6{x — x'{l — z)) 



a /•! , 1 rfr^°™(x') , X 

dx ; (1 

2iT Jx x' dx' x' 
-Fo C dx'{x' ~ x)('i-2x^\ 

Tl Jx ^ ' 



a 



-ro((l-a;)2(5-2a;)) , (23) 

where dV^'^'^x') j dx' = To2x'\3 - 2x'), and 6{x - x'{l - z)) = 6{z - ^)/x'. The 6- 
function expresses energy conservation in the splitting process. The integration over z 
shifts the lower boundary of the x' integration from to a; because of the (5-function . 
The final result exactly reproduces the result 

Numerically, the flip spectrum function is rather small compared to the 0{a) no- flip 
spectrum function. However, when averaging over the spectrum the ratio of the 0{a) flip 
and no-flip contributions amounts to a non-negligible (—12%), due to cancellation effects 
in the 0{a) no- flip contribution. It is clear that the corresponding helicity flip effect is 
larger in QCD due to ag/a ~ 10. 
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